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Adsorption of Gases on Heterogeneous Solid Surfaces. A New Method for Evaluating 
the Differential Distribution of Adsorption Energy, Based on the Polynome 
Approximation of the Distribution Function 


WSTĘP 


Zagadnienie heterogeniczności adsorbentów jest interesujące zarówno 
z teoretycznego, jak i praktycznego punktu widzenia. Problemowi temu 
zostało poświęconych wiele prac z dziedziny adsorpcji fizycznej [1—5], 
lecz istniejące rozwiązania ciągle cechuje pierwsze, niezbyt dokładne sta- 
dium. Pierwszą trudnością, jaką napotykamy przy opisie adsorpcji na po- 
wierzchniach heterogenicznych jest problem ilościowego opisu hetero- 
geniczności adsorbentów. Jeśli założymy, że struktura adsorbentu nie ule- 
ga zmianom wraz ze zmianami ciśnienia adsorbatu, wtedy pełnym iloś- 
ciowym opisem własności energetycznych adsorbentu jest rozkład energii 
adsorpcji e w przestrzeni adsorpcyjnej. Innymi słowy, należy określić 
funkcję [«(r), gdzie r jest wektorem położenia w przestrzeni adsorpcyj- 
nej lub ogólnie — zespołem zmiennych w statystycznej przestrzeni fazo- 
wej. Określenie funkcji [e(r)] z danych eksperymentalnych jest możliwe 
tylko w przypadku powierzchni homogenicznej. W przypadku powierzchni 
heterogenicznych ogólnie przyjętą ilościową charakterystyką własności 
energetycznych adsorbentu jest funkcja zdefiniowana następująco [3]: 
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av 

x (e) = rr (1) 
gdzie V jest objętością przestrzeni adsorpcyjnej lub ilością centrów 
adsorpcyjnych, w zależności od przyjętego modelu adsorpcji (adsorpcja 
niezlokalizowana lub zlokalizowana). Funkcja [z(s)] nazywana jest róznicz- 
kowym rozkładem energii adsorpcji. Proponowano wiele metod wyzna- 
czania różniczkowego rozkładu energi adsorpcji z danych eksperymental- 
nych. Pomijając metody kalorymetryczne i spektroskopowe, pozostałe 
mogą być podzielone następująco: 

1) metody wykorzystujące zależność temperaturową danych adsorp- 
cyjnych; 

2) metody wykorzystujące zależność ciśnieniową danych adsorpcyj- 
nych. 

Pierwsza grupa metod opiera się na wirialnym formalizmie opisu zja- 
wiska adsorpcji [5, 6]. Druga grupa metod bazuje na założeniu płatowego 
modelu adsorpcji. Całkowitą izotermę adsorpcji ©, opisuje sie tu nastę- 
pującym równaniem całkowym: 


6: (p, T)— J 0. (p, T, e) x (e) de (2) 
2 


gdzie ©, jest lokalną izoterma adsorpcji, opisującą równowagę adsorpcyj- 
ną na homogenicznym płacie o energii e, zaś 2 jest przedziałem, w któ- 
rym zawierają się energie adsorpcji istniejące realnie na powierzchni. Za- 
gadnienie wyznaczania funkcji ze sprowadza się do rozwiązania równa- 
nia całkowego (2). Spośród trzech występujących w równaniu (2) funkcji 
z danych eksperymentalnych określa się tylko funkcję O. Podstawowym 
założeniem czynionym w celu rozwiązania równania (2) jest przyjęcie 
określonego modelu adsorpcji na homogenicznym płacie powierzchni, in- 
nymi słowy, określenie postaci analitycznej funkcji Oi. 

W celu wyznaczenia funkcji z(e) należy rozwiązać równanie (2). Istnie- 
ją dwie alternatywne drogi postępowania: 

1) założenie a priori określonej postaci funkcji (e), zależnej od para- 
metrów, obliczenie całkowitej izotermy adsorpcji i po zastosowaniu pro- 
cedury best-fit określenie liczbowych wartości parametrów w funkcji roz- 
kładu [2, 8]; 

2) analityczne rozwiązanie równania (2) po uprzednim założeniu pew- 
nej postaci analitycznej dla 0. 

Osobną grupę metod stanowią graficzne metody Hobsona [9, 10] 
oraz Adamsona i współprac. [11, 12]. 

Celem niniejszej pracy jest zaproponowanie pewnej nowej metody 
wyznaczania funkcji z(e) z danych wagowych adsorpcji. 
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Wzorem wielu autorów, przyjmować będziemy, że równowaga adsorp- 
cyjna na homogenicznych płatach powierzchni opisana jest równaniem 
Langmuira. Przyjmijmy więc, że 9, dana jest równaniem: 


1 
1 + (K/p) exp (—e/RT) 


gdzie K jest stałą związaną z funkcją sumy stanów cząstki zaadsorbowa- 
nej, ze względu na wibracje, rotacje i wewnętrzne stopnie swobody mo- 
lekuły zaadsorbowanej w warstwie pierwszej [13, 14]. Z prac Kiseleva 
[15] i Hobsona [9] wynika, że z dobrym przybliżeniem stałą K można 
wyrazić równaniem: 


9, (p, T. e) (3) 


K=1,76 - (MT)? - 104 (4) 
gdzie M jest masa czasteczkowa adsorbatu. Równanie (2) przybierze za- 
tem postaé: : 


X (e) de 
1 + (K/p) exp(—e/RT) 


9,9. = | (5) 
9 

Najeze$ciej spotykanym sposobem rozwiazania równania calkowego (5) 

jest metoda polegajaca na doprowadzeniu tego równania do postaci tran- 

sformaty Stieltjesa [16—20]. Wprowadzajac bowiem notacje: 
K/p—y—1; exp (e/RT)=x+1 (6) 

i zakładając, iz obszar zmienności e, Q — «C0, oo, otrzymamy wówczas tran- 

sformate Stieltjesa w postaci: | 


0, (; K | = aor + VI, 


— ,T P (7) 


1 


0 
Stosując równanie (7) Sips [16, 17] udowodnił, że rozkład eksponencjal- 
ny: 


` x(e)=[exp (6/RT)—1]-© O<C<1 (8) 
odpowiada uogólnionej izotermie Freudlicha: 
\c 
0, (p, T)= (DK (9) 


Najcenniejszy rezultat na bazie tej metody został otrzymany przez M i- 
srę [18—20], który wykazał, że izoterma globalna Dubinina—Radusz- 
kiewicza: 

0.(p,T)=exp [—B {RT In(K/p+1))?] (10) 


która, jak stwierdzono, dobrze opisuje adsorpcje na wielu heterogenicz- 
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nych adsorbentach [21—22] koresponduje z prawostronnie rozmytym roz- 
kladem gaussowskim: 
y(e) =2 Be exp ( Be?) (11) 
W niniejszej pracy problem wyznaczania funkcji dystrybucji energii 
adsorpcji rozwiazano w sposób nastepujacy. Rozwijajac izoterme Lang- 
muira, dana równaniem (3) w szereg Taylora wobec p otrzymujemy: 


i SAN A 25 
RATA * e (12) 


Korzystając z powyższego przedstawienia, równanie (5) możemy zapisać 
w następującej postaci: 


6,(p,T)= Ya WE | "exp T | x (e) de (13) 
2 


Funkcja dystrybucji energii adsorpcji z(s) może być z powodzeniem 
aproksymowana wielomianem l-tego stopnia: 


L 
1 (e = ŻA, gi (14) 


Podstawiając rozwinięcie (14) do równania (13) otrzymujemy: 


— l 
8,(p, T) Di () (x) | exp lex DZE de (15) 
nl Q 


Problem znalezienia funkcji z(e) sprowadza się zatem do wyznaczenia 
współczynników 4, wielomianu danego równaniem (15). W tym celu zróz- 
niczkujemy n krotnie obie strony równania (15) względem p. Biorąc na- 
stępnie wartości pochodnych w punkcie p=0, otrzymujemy: 


l 
9^6, EP UNE. ne i 
paka sl 1) Kn Ve (75) M At (16) 
2 i=0 


Przedział możliwych zmian energii O jest zazwyczaj określany a prio- 
ri. Wielu autorów zakłada przedział zmienności = , oo) [2, 8—10, 
16—20], choć faktycznie energie adsorpcji zmieniają się w pewnych okreś- 
lonych granicach od minimalnej wartości emin do maksymalnej wartości 
Emax energii [2, 3, 23, 24]. W dalszych naszych rozważaniach przyjmować 
będziemy, że obszarem możliwych zmian e jest przedział domknięty 
Q=<0, A>. W celu ułatwienia zapisu wprowadźmy następujące ozna- 
czenia: 


A 
Fr er. KV. i (17) 
0 
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mo 
B 1 ( t 


m : 
m m! << F s GR 


Z równań (16), (17), (18) otrzymamy wówczas: 
Bm = AF. o AyF mat ... + AF na (19) 


Wyznaczenie współczynników A, jest już łatwe. Wystarczy w tym celu 
rozwiązać następujący układ równań liniowych: 


Bn=ApF mot AiF, i ... AF mi 

Bm+1= AqgFm+1,0 T A1 Fi. i APR 

B,— ApF pot AI Fo.i ... AF 

Buica AFN o AF miga... tA instar 


Współczynniki B, mogą być wyznaczone z danych eksperymentalnych 
przez kolejne graficzne bądź numeryczne różniczkowanie początkowego 
odcinka izotermy, zaś współczynniki Fp obliczamy, traktując początko- 
wo A jako parametr. Liczbową wartość A obliczamy przez zastosowanie 
procedury best-fit z następującego równania: 


Baa AF AiF ... + AF mitt (21) 


posługując sie współczynnikami A; obliczonymi z równania (20). Przez 
porównanie wielkości liczbowej Bus obliczonej z izotermy doświadczalnej 
i wyznaczonej z równania (21) korygujemy wartość liczbową współczyn- 
nika A aż do uzyskania najlepszej zgodności danych eksperymentalnych 
i teoretycznych. W ten sposób po dokonaniu szeregu powyższych ope- 
racji uzyskujemy numeryczne wartości współczynników wielomianu (14). 
Współczynniki te możemy określić następującymi równaniami: 


(20) 


| B. A Sa Fn+11 

Bn+1 Fm 1. 1 — m+l—1,1 

B, Epi LE A (22) 

AD= Barisi Fm 1-1. i gi" im m Fm+t-1, 
D M 
jg * F Ba Fmi+1 e tors Fn 
Fm+1,0 ELIT Fma3,—1 Bn+1 F mitta ale be Fn+1a 
Fy,o Fp. 1 B, Pos Foz (23) 
AN = Fm+i-10 — Fmns- Pe Fn+t1,1+1 · Li. Faria 


8 — Annales 
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gdzie 
Em 0 Fo Fh 
ER E431 Fanti 
M=|F,, By aU Rl (24) 
rae Facta t Frisia 


Wskazniki j oznaczaja, ze wspólczynniki A, obliczone zostaly po uprzed- 
nim dopasowaniu parametru. 

Zgodnie z powyzszymi równaniami funkcje dystrybucii energii ad- 
sorpcji opisuje następujące wyrażenie: 


i 
X(e) = Y Ape (25) 
i=0 


WYNIKI NUMERYCZNE I DYSKUSJA 


Do przeprowadzenia obliczeń ilustracyjnych użyto dobrze znanej w li- 
teraturze izotermy adsorpcji argonu na rutylu, zmierzonej przez Draina 
iMorrisona [25]. Obliczenia przeprowadzono przy użyciu funkcji dy- 
strybucji energii adsorpcji danej równaniem (14), kładąc 1=2. 

(e) = A Ayet Age? (26) 

Wyznaczenie funkcji dystrybucji danej równaniem (26) sprowadza się za- 
tem do rozwiązania następującego układu równań liniowych: 

BI = AF + AF, ¡+ A, 

B= AF2, + AF21 + A2F. 2 (27) 

B = A F3,9 + AyF 3,1 + AoF 3,2 
Współczynniki Bi, Bs, Ba, By, obliczono przez graficzne rózniczkowanie 
doświadczalnej izotermy adsorpcji. Wynosza one kolejno: B,=1,026 
25 108, B = 1,248 03 - 1012, B,— 1,613 68 - 1015,B,— 1,85 330 - 10%, zaś współ- 
czynniki F sa określone w tym wypadku następująco (p=1, 2, 3): 


A 
F7 C (5) de T exp (Ar- | (28) 
0 7 s 


å 
po żel eb (PR -2T | (24). TRU (et) RT| (29) 
pA | p 25) sas p exp 4 mo RT)" p 
0 
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_ RT [24 alma 17. _RT | 2(RTY 
-PT exp(P4). 42 =. exp [22 A : exp (z TUR p° (30) 


Parametr 4=2500 cal, obliczono poprzez zastosowanie procedury best-fit, 
sposobem podanym w części teoretycznej z równania: 

B. AFA AFA AF. (31) 

Uzyskano w ten sposób nastepujaca funkcje dystrybucji energii adsorpcji: 

xke) = — 7,833 52-10-8 e 3, 211 74- 10-4 - 3,1334 - 10-1 (32) 


Xt£) 
(cal) 


1500 2000 2500 3000 3500 S (cal) 
Ryc. 1. Funkcje dystrybucji energii adsorpcji dla ukladu argon na rutylu [25]; 
1— dystrybucja DR wyliczona z równania (11) B=9,829 10—7 (cal/mol)-?, 2— dy- 
strybucja wyliczona z równania (32), 3— dystrybucja wyznaczona przez Adam- 
sona [26] 


Obliczona przez nas funkcje dystrybucji energii adsorpcji przedstawia 
ryc. l. Funkcja ta zostala porównana z rozkladem DR, danym równa- 
niem (11), i z rozkładem Adamsona [26]. Mimo użycia tylko drugiego 
stopnia wielomianu (26), zgodność obu funkcji jest zadowalająca. 

Zalety proponowanej tu metody w porównaniu z innymi są następu- 
jące. Po pierwsze nie występuje konieczność zakładania określonej anali- 
tycznej postaci funkcji O., bowiem do wyznaczenia funkcji z(e) wystarcza- 
ją w pełni dane eksperymentalne. Po drugie funkcja dystrybucji energii 
adsorpcji wyznaczana jest z początkowego odcinka izotermy eksperymen- 
talnej, kiedy założenie związane z zaniedbaniem wzajemnych oddziały- 
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wan molekuł adsorbatu (model Langmuira), jest najlepiej spełnione. Sama 
zaś funkcję %% można przybliżać wielomianem w sposób dowolnie do- 
kładny. 
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PE3IOME 


B nammoú pa6ore npencragnen HOBBI Merog onpenenenun nucpchepen- 
IMaJIbHOTO pacIipeqeJleHMA SHepruu axcopOnum. JTOT Mero] OCHOBaH Hà 
JIOCKYTHOŃ MOJĄEJIM HEOAHOPOAHOÓ HIOBepXHOCTM M Hà MHOTOYJIEHHOŃ an- 
HpOKCHMAHNUM dyHkKUMM paciipereJleHMA. ; 
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SUMMARY 


A new method of evaluating differential distribution of the adsorption 
energy is proposed. The method is based on the assumption of the patch 
model of the heterogeneous surface and on the polynome approximation 
of the distribution function. 
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